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ABSTRACT 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



The modern maritime battlefield is dominated by the new 
generation of sea-skimming, high-speed, stealthy and highly 
agile anti-ship missiles. Anti-ship cruise missile 
technology continues to evolve and to overcome the 
performance of the existing ship self-defense weapon 
systems. The Free Electron Laser (FEL) could be the 
ultimate speed-of-light , hard-kill weapon system, offering 
unique features such as tunability, high power, pinpoint 
accuracy and infinite magazine. Multimode computer 
simulations were used to explore the operation of the Thomas 
Jefferson National Acceleration Facility (TJNAF) FEL with 
untapered and positively tapered undulator. The final 
steady-state power, the steady-state gain and the electron 
energy spread as a function of desynchronism were determined 
for both 34.5 Mev and 47.5 Mev electron beam energies. 

This thesis also includes an experimental study of 
damage induced to Polyimide Fiberglass and F2 Epoxy samples, 
by the TJNAF FEL. Irradiations of the samples were 
conducted changing various parameters, such as the 
wavelength, average power, pulse repetition frequency, cross 
wind and spot size in order to explore the damage mechanism. 
At this stage of evolution, TJNAF FEL is now capable of 500W 
output average power. In order to achieve the required 
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to a small 



2 

intensity of 10 kW / cm the beam was focused 
radius, and scaling guidelines were developed 
predict the damage caused by a high power laser 
area . 



in order to 
over a large 
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EXECUTIVE SUMMARY 



The Free Electron Laser (FEL) could be the solution to 
the complex problem of ship's self defense against modern 
anti-ship missiles (ASM). The FEL ' s immediate response, 
infinite magazine, speed-of-light target engagement, and 
continuous tunability make it an excellent alternative of 
the current Close In Weapon Systems (CIWS) . 

Computer simulations were used to explore the operation 
of the Thomas Jefferson National Acceleration Facility 
(TJNAF FEL) with untapered and positively tapered 
undulators . Taper rates of 5=+4n, +6n, and +8n were used. 
The final steady-state power, the steady-state gain in weak 
fields, and the electron energy spread as a function of 
desynchronism were determined for both 34.5 Mev and 47.5 MeV 
electron beam energy. The shape of those curves provides 
useful information about the FEL operation. The final 
power and gain were significantly smaller for the 47.5 MeV 
electron beam than for the 34.5 MeV electron beam, because 
of the lower value of dimensionless current j . However, the 
general shape of the curves is very similar for each energy. 
The desynchronism curves present a sharp peak at small 
values of desynchronism accompanied by the trapped-particle 
instability. The gain versus desynchronism curves peak at 
moderate values of d and give near few gain for both small 
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and large values of desynchronism. Limit cycle behavior was 
observed in the final power for both the 34.5 MeV and 47.5 
MeV electron beams, but not for all tapers. The steady- 
state power and the steady-state gain turned out to be 
higher for the conventional undulator than for the tapered 
undulator. Furthermore, the tapered undulator induces more 
energy spread to the electron beam, and demonstrates less 
efficiency for beam recirculation. The width of the 
desynchronism curve decreased as the tapering rate is 
increased. The operating width of the desynchronism curves, 
Ad, as a function of taper has been measured in experiments 
conducted with the IR Demo FEL at Jefferson Lab [Ref 16] . 
Simulation results are summarized in Table 1. There is good 
agreement between theory and experiment. 



Electron Beam 

Energy 


34.5 MeV 


47.5 MeV 


Taper- 6 


0 


4n 


6n 


8n 


0 


4n 


6n 


8n 


Maximum Power 

(dimensionless) 


206 


141 


96 


98 


170 


112 


89 


9 


Maximum 
Gain ( % ) 


85 


58 


33 


12 


76 


53 


32 


12 


Desynchronism 

Curve 

Width ( 

x NX) 


0.38 


0.3 


0.19 


0.06 


0.34 


0.26 


0.18 


0.04 


Limit Cycle 
Behavior 


YES 


NO 


NO 


YES 


YES 


YES 


NO 


YES 


Maximum 

Electron Energy 
Spread (%) 


8 


8 


7.5 


7 


7 


6.5 


6 


2 



Table 1. Summarized Simulation Results. 
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This thesis also includes an experimental study of 
damage induced to missile nose-cone materials by the FEL. 
Experiments were conducted with the TJNAF FEL in order to 
study the damage caused by the laser beam to some sample 
materials commonly used in missile cones. The TJNAF FEL is 
capable at this point of 500 W output power. Since 1 MW 
power over a spot of 100 cm 2 is approximately required to 
shoot down a missile, the beam was focused to a small spot 
area in order to achieve the required intensity of 10 

kW / cm . In this way several guidelines were developed that 
could predict the damage caused by a high power laser over a 
larger area. The pulse repetition frequency (PRF), the 
wavelength, the average power, the spot size and the 
crosswind speed were varied during the irradiations in order 
to determine the influence of these parameters on the damage 
induced. Airflow of 60-80 mph speed did not have a 
significant effect, and increased the penetration depth rate 
by 10%-15%. The penetration depth rate is proportional to 
the laser intensity. The A=3.1 ym wavelength was more 
effective than the A=4.875 ym wavelength. In addition, the 
lower PRF (18.7 MHz) seems to contribute to higher 
penetration rates due to the higher energy per pulse. 
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I. 



INTRODUCTION 



In recent years highly sophisticated anti-ship missiles 
(ASMs) have been introduced to the sea warfare environment 
making the modern maritime battlefield much more compressed 
and deadlier than ever before. The ability to react quickly 
and correctly to multiple ASM attacks is now an essential 
part of the self-defense capability of a warship. Current 
ship weapon systems tend to lose their ability to meet this 
requirement as the offensive missile technology evolves 
spectacularly. A Free Electron Laser (FEL), as a speed of 
light, hard-kill weapon system, is a potential solution 
offering many attractive characteristics in performance, 
readiness, and life cycle management. 

Chapter II of this thesis describes the FEL, and its 
principles of operation, and gives an overview of background 
theory. Chapter III addresses the benefits of the FEL as a 
hard kill weapon system against anti-ship missiles, and 
analyzes the shortcomings and needs of current ship self 
defense . 

Chapter IV describes the operation of the tapered 
undulator design, and simulates the performance of the 
Thomas Jefferson National Acceleration Facility (TJNAF) FEL. 
Multimode simulations are used to analyze the operation with 
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untapered and positively tapered undulator, describing the 
evolution of short optical pulses at far infrared 
wavelengths. Trapped particle instability and limit cycle 
behavior are observed and analyzed. Furthermore the steady 
state power, the steady state gain and the electron energy 
spread as a function of desynchronism are determined for 
both 34.5 Mev and 47.5 Mev electron beam energies. 

Chapter V presents the analysis of experimental 
irradiations on several materials commonly used in missile 
bodies. The goal of the experiments was to determine the 
amount of damage induced by the unique FEL beam of short 
pulses and develop scaling laws that could predict the 
damage caused by a much higher laser power over a large 
area . 
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II . 



FREE ELECTRON LASER THEORY 



A. DESCRIPTION AND BASIC OPERATION 

The Free Electron Laser (FEL) concept was first 
introduced by John Madey in 1971. The basic parts of a FEL 
are the electron accelerator, the undulator, and the optical 
resonator cavity as shown in Figure 1. 




The undulator is a configuration of magnets or 
superconducting coils, which produce a spatially periodic 
magnetic field. A relativistic electron beam generated by 
the accelerator is injected into the undulator. As 
electrons travel through the undulator, they are 
periodically deflected in the transverse direction by the 
imposed magnetic field. This transverse motion causes the 
electrons to radiate (spontaneous emission) in a narrow 
forward cone. The electromagnetic radiation is stored in 
the optical cavity and the optical power grows over many 
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passes inside the resonator. The coupling of the optical 
field, the magnetic field, and the electron beam allows 
extraction of kinetic energy from the relativistic electrons 
and conversion to high optical power. Subsequent electrons 
moving inside the undulator interact with the continuously 
growing optical field and radiate in the presence of light; 
this is stimulated emission. 

The combined effect of the undulator magnetic field and 
the electromagnetic radiation gives rise to the axial 
bunching of the electrons at the optical wavelength, 
allowing them to radiate coherently. 

B. PARAMETERS AND TYPICAL VALUES 

The FEL ' s most important operation parameters are 
listed in Table 2. 



PARAMETER 


DEFINITION 


RANGE 


B 


Undulator Magnetic Field (rms) 


2-7 KG 


N 


Number of Undulator periods 


- 100 


^0 


Undulator Wavelength 


2-10 cm 


X 


Optical Wavelength 


5nm-5cm 


L 


Undulator Length 


1-20 m 


K 


Undulator Parameter 


0.1-10 


Y 


Lorentz Factor 


5-5000 


ymc 2 


Electron Beam Energy 


MeV-GeV 


i 


Electron Beam Current 


1-1000A 


r b 


Electron Beam Radius 


1mm- 1cm 


Pin 


Input Power 


10-100 MW 


P out 


Output Power 


1-10 MW 


n 


Efficiency of Energy transfer 


l%-50% 



Table 2. FEL's Operation Parameters. 
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The undulator parameter K is defined by K = eBA 0 / 2nmc 2 
and its value is usually around 1. The radiation optical 
wavelength is A = A 0 (l + K 2 )/2y 2 (Ref. 1) . The FEL ' s 

operating wavelength can be easily tuned simply by adjusting 
the electron beam energy or the undulator parameter. It is 
worth noting that the FEL can be designed to work in an 
extremely wide spectrum of wavelengths and in new wavelength 
ranges, where there are few powerful sources of radiation 
like 20ym to 500ym and 20nm to 300nm. This remarkable 
feature makes the FEL an attractive device, demonstrating 
its potential for various applications. 



C. ELECTRON EVOLUTION INSIDE THE UNDULATOR 

1 . The Undulator 

Two kinds of undulator polarizations are commonly used 
in FELs, the helical and the planar. The helical undulator 
consists of superconducting coils and generates a magnetic 
field B = B [cos(k 0 z) , sin(k 0 z) , 0] on axis with k 0 = 2n / A 0 , 
and it produces circularly polarized light. The planar 
undulator consists of an arrangement of magnets and 
generates a magnetic field B = B [o , sin(k 0 z) , o] on axis, 
producing linearly polarized light. 
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Figure 2. Helical Undulator. Figure 3. Linear Undulator. 

(From Ref. 2) . (From Ref. 2) . 



2 . Resonance Condition 

As electrons move through the undulator in the axial 
direction, they are periodically deflected in the transverse 
direction by the spatially varying magnetic field. A 
typical motion of electrons inside the undulator is shown in 
the figure below. Note that z and x are the axial and 
transverse directions respectively. 




Figure 4. Electron motion inside the Undulator. 
( From Ref . 2 ) . 
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As an electron moves through an undulator wavelength 
A 0 , it completes a full oscillation period in the transverse 
direction and emits a wavelength of light A. So one 
wavelength of light traveling with speed c passes over a 
slower electron, as it travels through an undulator 
wavelength A 0 , with speed (3 z c . This particular electron- 
photon "race" is called the resonance condition and is shown 
schematically in the figure below. 




The electron acceleration causes the emission of 
synchrotron radiation in a cone of width 1/y (Ref. 2), 

where y = 1 / -/l - (J 2 is the Lorentz Factor (3 = u/c, and u is 
the electron's speed. 
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We define (3 2 c the electron speed in the axial 



direction, where c is the speed of light. We can easily 
calculate the wavelength A of the emitted radiation using 
the electron-photon speed difference. 

( \ 



A 0 1 


Xp, / \ 


f 1 ) 




= -2. • (l - pj = A 0 


1 


TO 

N 

o 


Pz 





2y 2 J 



*0 



2y 2 - 1 




and 



* P. 



For highly relativistic electrons, y >> 1 and 

A * A 0 / 2y 2 . 

3 . Pendulum Equation 

We consider a helical undulator magnetic field 
B = B [cos(k 0 z) , sin(k 0 z) , 0], where k 0 = 2n/A 0 is the 

undulator wave number and A 0 is the undulator wavelength. 
We assume for the optical field a circularly polarized plane 
wave, with electric and magnetic fields 

E s = E [cos'? , -sin'F , o] and B s = E [sin 1 ? , cos'? , o] respectively 
where W = kz-cot + cp , <p is the optical phase, 

k = 2n/A is the optical wave number, and A is the optical 
wavelength. It should be pointed out that the electric and 
magnetic fields E s , B s have the same amplitude E, as they are 
measured in cgs units. 
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Relativistic electrons with charge magnitude e are 
injected in the undulator with energy ymc 2 and momentum 
P = ymu , where y is the Lorentz factor, m is the electron 



mass, u = Pc = p ± c + P z c is the electron's velocity, and 



Pj_c, P z c are the velocity components in the transverse and 
axial direction respectively. Electrons interact with the 
fields and they experience relativistic force 



dP 

F = — = -e 
dt 



From [ 1 ] 



_ u - 
E + — x B t 
c 



d(yp) 

dt 



[e + P x B t ], 



:i] 



[ 2 ] 



me 



where B T is the total magnetic field of the undulator and 
optical field combined together § T = B + B s . The magnetic 
field is always perpendicular to the electron motion and 
does not affect the electron ' s energy ymc 2 . The electron 
energy change is given by 



d(ymc ) - - _ dy e - - 

= F • u = — eE • u => — — p • E . 

dt dt me 



: 3 ] 



Inserting the fields into the Lorentz equation we get 



= — [E ( 1- p z ) • (cos % 
dt me 

^ [0, 0, Efccosf 
dt me 



- sin T, 0 ) + p z B • (- sin (k 0 z) , cos (k 0 z) , 0)] 

- P y sinT) + B • (p x sin(k 0 z) - P y cos(k 0 z})] 



[4] 

[5] 
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— = - — E[p x cos T 


- p y sin w]. 




[6] 


dt me L 






For relativistic 


electrons 


E(1 - P z ) « Bp 2 , 


so we can 


disregard the term 


E(1 — p z ) • 


(cosT , - sinT , 0) . 


Within this 



approximation, integration of equation [4] with respect to 
time gives the transverse motion as 



— K 

Pi [cos (k c z) , sin ( k c z) , o] 

Y 

where K is the undulator parameter and 



(3 X cos (k 0 z) , 

Y 



K 



p = sin (k 0 z) . 



[7] 

[ 8 ] 
[9] 



We substitute P x , p y into equation [6] and obtain 
dy . eKE . v 

— = Y = cos(^ + <p). [10] 

at ymc 

The Lorentz Factor y describes the electron's energy, 
£ = (k + k c ) z - at is the electron's phase generated by the 
combination of the undulator and optical fields, and <j> is 
the optical phase. Since A « A 0 , k » k 0 and y»l 

(relativistic electrons) , the initial electron phase in the 
beginning of the undulator at t=0 is £ 0 = (k + k 0 )z 0 = kz 0 . 
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Combining this with 



From equation [7], we have = — . 

Y 2 



the Lorentz factor, we get 



f 2 = 1 - = 1 - g -pi-1 *l-g. 

r 



K z 



Differentiation of [11] with respect 

i = 

y (l + K 2 ) 



[ 11 ] 

to time gives 

[ 12 ] 



which relates y to P 2 . Differentiating ( twice, we acquire 



4 = (k + k 0 )z - to = (k + k 0 )cP 2 - (o 



and l = (k + k 0 )cP 2 => P 2 = 






[i3; 



[14] 



(k + k 0 )c 

We substitute equation [14] into [12] and get 






[15] 



y (l + K 2 ) • (k + k 0 )c ' 

We now make use of the resonance condition 
A = A 0 (l + K 2 )/ 2y 2 or k 0 = k(l + K 2 ) / 2y 2 , and assume that for 
relativistic electrons P 2 = 1 and k + k 0 = k, since k » k 0 . 
Equation [15] becomes 






y 2k 0 c 



[16] 



Combining equations [12] and [16], and solving for 4 we 
obtain 
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[ 17 ] 



2k 0 eKE , v 

C = — i — cos (< + ‘P) ' 

y m 

which describes the dynamics of the electron phase by the 

Pendulum Equation . 

The time interval for an average electron to cross the 
undulator is At s L / |3 0 c , which motivates the introduction 

of the relevant dimensionless time t = 3 0 ct / L = ct / L , 
ranging from 0 to 1 along the whole undulator length L. We 



OO 

use the notation of £ as the second derivative of £ with 

respect to x, and we define |a| = 4nNeKLE / y 2 mc 2 as the 

dimensionless optical field amplitude. The pendulum 

equation takes its final form 

OO 

£ = M cos ( £ + cp) . [18] 

Certain important remarks must be made from this analysis: 



• £ is a 


microscopic 


variable describing 


the 


relative 


position of 


an 


electron within 


an 


optical 


wavelength, as 


it 


travels through 


the 



undulator. 

• Equation [10] relates y , which is the electron's 
rate of energy change, to the electron's phase 
(^ + cp). As cos(( + cp) takes values from -1 to 1 
within one optical wavelength, half of electron's 
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phases result in y > 0, and the other half give 
y < 0, which implies that electrons become spread 
in energy. 

• The fact that y cc EK / y implies that stronger 
fields result in a larger electron energy change, 
whereas more energetic electrons (larger y) 
reduce the coupling of the undulator and optical 
fields . 

• Electrons becoming spread in energy within an 
optical wavelength, implies that some of them 
gain energy and move faster while other lose 
energy and move slower, thus they become bunched 
on the scale of the optical wavelength. 

• Electron motion inside the undulator, which is 
rather complicated, is described by the well- 
known pendulum equation, which facilitates 
remarkably the visualization of the electron 
phase dynamics, and the understanding of the 
phase space diagrams. 
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D . DIMENSIONLESS PARAMETERS 

Several dimensionless variables can summarize recurring 
combinations of physical parameters, which govern the FEL 
operation, and are used to gain insight into the relevant 
physical processes. An important quantity characterizing 
many FEL properties is the undulator parameter 

K = eBA 0 / 2nmc 2 where B is the rms field strength over each 
period, A 0 is the undulator wavelength, and e, m, are the 

electron charge and mass respectively. In most FELs K ^ 1 . 
A pass of the electron pulse through the undulator is 
described by the dimensionless time x = ct / L where 
t = 0 * 1 . The electron phase £ and phase velocity 

v = ^ = L[(k + k 0 ) • (B z — k], describe the electron velocity and 

position in a section of the beam one wavelength long, and 
define bunching and bunching rate of the electrons. The 
phase and phase velocity are coupled by the dimensionless 

complex optical field a = |a|e 1<p through the simple pendulum 

OO 

equation £ = |a|cos( ( + cp) . The field amplitude is 

|a| = 4nNeKLE / y 2 mc 2 , and determines the electron-bunching 
rate on the optical wavelength scale. When a « n the 
optical field is weak, and when a » n the FEL operates in 
the strong field regime. The response of 
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the optical wave to bunching in the beam is described 
by the dimensionless current j = 8N(enKL) 2 n e / y 3 mc 2 , where n e 

is the electron density within the electron beam. The 
dimensionless current provides the coupling between the 
electron beam and the light wave through the wave equation 

0 7 

a = -j(e y (Ref. 1) . When j«n the FEL gain is low, and 
when j»n the FEL gain is high. Each definition above, £ 
and v describing the electrons, a describing the light 
wave, and j describing the whole FEL, is physically 
meaningful and useful in evaluation of many FEL effects. 

E . GAIN SPECTUM 

Figure 6 shows the FEL gain spectrum for a typical low 
current FEL with j=2 and N=21, as the initial optical field 
amplitude a 0 evolves from 0 to 40. G is the FEL gain for 

one pass of the electron beam through the undulator in a 

single optical mode, and v 0 is the initial electron phase 
velocity at the beginning of the undulator. Figures 7, 8 

and 9 present a cross section of the gain spectrum for 
a 0 = 0, a 0 = 20 and a 0 = 40 respectively. The plots were 
generated by computer simulations. The gain has a maximum 
value (peak gain), for each optical field strength a 0 . 

Furthermore, the position of the peak gain changes as a 0 

increases, and the gain spectrum distorts. 



15 




Figure 6. FEL Gain Spectrum, j=2. 
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The value and the position of the peak gain are crucial 
to the FEL operation, since they determine the optical 
frequencies that are going to be amplified as well as the 
final optical power. The FEL gain is clearly antisymetric 

about resonance v o =0. In weak optical fields the gain 
spectrum has a peak value of G « 0.13j at v 0 = 2 . 6 , and 
width Av o « n . The gain is zero at resonance. As the 

optical field grows to saturation the gain spectrum 
distorts, and becomes broader with width Av 0 « Va . The peak 
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gain decreases substantially and moves to higher values of 
v 0 . In the case of j = 2, the peak gain drops from 27% in 
weak optical fields, to 1% when the optical field saturates. 
The phase velocity of peak gain moves from v 0 =2.6 to 
v o =10.8. The corresponding change in the wavelength is 
AX/A = Av 0 /4nN = 3% in this example. 

F . FEL ' S APPLICATIONS 

The FEL mechanism can be designed to operate over a 
large region of wavelengths, from nanometers to centimeters. 
Each FEL also has the capability of being continuously tuned 
simply by changing the electron beam energy. Most often the 
FEL is operated so that the net energy transfer is from the 
electron beam to the radiation field in order to amplify it. 
Nevertheless, it is possible to operate the FEL so that the 
net energy transfer is from the radiation field to the 
electrons making an accelerator or an inverse FEL. 

The FEL has become an exciting alternative to other 
radiation sources, like microwave tubes and lasers, and can 
extend the operational range of both. Because of its 
flexibility, the FEL can be applied in many diverse areas 
such as lithography, plasma heating, particle acceleration, 
as well as material, biological, medical, and solid state 
research . 
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III. THE FEL AS SHIP SELF- DEFENSE WEAPON 



A. THE ANTI -SHIP MISSILE THREAT 

In the modern naval battlefield, anti-ship missiles 
(ASM) have become a very serious threat even to the most 
modern and complex warships. Faster and faster ASMs have 
dramatically decreased response times for surface vessels. 
A simple scenario helps us realize the complexity of the 
problem. 

Consider a sea skimming ASM, flying less than 10 meters 
above sea level and heading towards a ship. Let the height 
of the ship's sensors be 25 meters above the sea surface, so 
that the horizon distance in which the missile can be 
detected and tracked is approximately 10 nautical miles 
(nm) . Given the missile's speed to be 1 mach = 600 knots = 
10 nm/min, it takes one minute for the missile to impact the 
ship. The time is decreased to 30 seconds if the missile 
approaches with speed mach 2. 

The Aegean Sea, where the Hellenic naval combatants 
usually operate, has a very specific geography. It is a 
complex of thousands of neighboring islands, less than 12 
nautical miles apart each other. This complicated battle 
arena gives rise to certain specific scenarios. For 
instance, a battleship operating 5 nautical miles from the 
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coast of a small island could suddenly be face to face with 
a Harpoon ASM emerging behind the island, using the waypoint 
navigation feature. 

Due to the spectacular development of anti-ship missile 
technology, ASMs have become state-of-the art weapons with 
highly supersonic speeds, ten's of g's maneuverability, 
multiple guidance modes, countermeasure capability, reduced 
radar cross section and IR signature, and pinpoint accuracy. 
There seem to be no end to this development making it hard 
for the ship self-defense systems to keep pace. 



B . SHIP SELF-DEFENSE 

The existing shipboard anti-missile systems are divided 
into two major categories: passive and active. Passive 

methods include decoys, chaffs, stealth technology, as well 
as tactics and techniques, which alter the ship's IR or 
RADAR signature, and generally deceive and confuse the 
missile sensors. 

On the other hand, active systems imply strong 
interaction with the incoming missile target and they are 
further divided into two categories: soft kill and hard 

kill. Soft kill methods disable the threat missile by 
attacking its electronics or blinding its sensor, but 



20 



otherwise leaving the target body undamaged. Hard kill 
weapon systems actually destroy the target-missile achieving 
large-scale physical damage. We are going to be focused on 
the hard kill method, which is the most desirable. A 
scenario of layered ship self-defense, involving both soft 
kill and hard kill methods, is demonstrated in Figure 10. 
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C. AVAILABLE HARD KILL WEAPON SYSTEMS 

Several hard kill weapon systems are used to confront 
the complicated ASM threat. 

1 . Large Caliber Guns 

The 5/54 or 3/62 OTOMELARA guns have effective ranges 
of 10 and 5 nautical miles respectively, and are the 
standard Hellenic Navy guns that can be used as a medium 
range defense, but with extremely low efficiency. The fire 
control system, the shell's speed, and the relatively slow 
rate of fire are not capable of shooting down an approaching 
supersonic sea skimming ASM with high maneuverability. 
Moreover, the shell dispersion (increasing with large range) 
limits the accuracy considerably and lowers the probability 
kill to minimal values. 

2 . Long Range Surface-Air Missiles 

Another approach considered as a long-range defense is 
the use of antimissile missiles. Several surface-to-air 
missiles, such as the standard Navy missile (SM-1, SM-2) or 

the Nato Sea Sparrow (NSSM) have been developed, and are 
equipped with excellent flight characteristics, 
maneuverability, and effective range. 

However, in order to take advantage of the long 
effective range, a long detection range is also required. 
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Since the sea skimming ASM missile's flight altitude is 
about 8 m with a corresponding detection range at the 
horizon of 11 nautical miles, making the long-range missile 
advantage becomes useless. Moreover, physical dynamics 
require that defensive missiles should have 2 to 3 times the 
maneuverability of the incoming missile. Current anti-ship 
missiles have the ability to fly at supersonic speeds, 
pulling more than lOg in terminal maneuvers. A simple 
calculation makes the maneuvering requirements of the 
defensive missile unfeasible. 

3 . CIWS 

The innermost layer of defense against anti-ship 
missiles is provided by the Close In Weapon Systems (CIWS) . 
The inner layer defense systems could be gun based, such as 
Raytheon's 20mm Phalanx, FABA's 20mm Meroka, Signaal's 30mm 
Goalkeeper, Oerlicon - Royal Ordance 35mm Millenium, or 
OTOBreda's twin 25mm Myriad. Some Navies advocate short- 
range missile based CIWS's such as the US-German RAM or the 
French Mistral. The basic characteristics and features of 
the currently most used CIWS are discussed below. 
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a) The Phalanx Gun 

At this time, the Phalanx Gun is the standard CIWS 
for both the U.S and Hellenic Navies. It is the most 
extensively used CIWS in the world with more than 700 
Phalanx systems deployed by 21 Navies. Since 1979, when it 
was initially deployed, the weapon has evolved from Block 0 
to Block IB through system improvements, upgrades, and 
modifications. The system is a rapid fire, computer 
controlled radar and gun system that automatically carries 
out search, detection, target threat evaluation, tracking, 
and firing. It is equipped with an electronically 
controlled, pneumatically driven, 20 mm six-barrel gun 
capable of firing 4500 rounds per minute within a range of 2 
km. Although the weapon provides superior performance 
against anti-ship missiles, it has several deficiencies. It 
experiences random and variable bullet dispersion, which 
diminishes its ability to destroy its intended targets at a 
range where missile debris does not hit the ship. The high 
rate of fire causes overheating of the barrels, so that the 
firing must be limited to bursts of only few seconds. The 
gun's magazine, which is capable of holding up to 1000 
rounds, can be expended in just four to five engagements and 
requires up to 30 minutes to be reloaded. If the missile 
attacks are more frequent, it means that the ship is left 
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with no self-defense. Finally, computer simulations have 
shown that the probability of hitting an incoming missile is 
extremely low for ranges more than 250 m (Ref. 11) . In 
Figure 11 we can see the probability of a single phalanx gun 
penetrator hitting an incoming missile flying at altitude of 
7 meters above sea surface with speed Mach 2, versus the 
missile range from the ship. It is easy to see that the 
probability drops very rapidly as the range increases. 




Figure 11. Probability of a single round to hit the missile 
versus missile range. 
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Figure 12, shows how the number of hits acquired 
by a missile with approach speed Mach 2 and altitude 7 m 
above the sea level, varies with range from the ship. The 
dashed line is for a firing rate of 3000 rounds/min (Block 
1A model) , while the solid line is for 4500 rounds/min 
(Block IB model) . Since 5 to 6 hits are needed to destroy 
the missile a typical killing range seems to be around 150 m 
from the ship. The missile may be destroyed within a range 
of 100 to 200 m, but its fragments contain enough kinetic 
energy to reach the ship and cause serious damage. 




Figure 12. Number of shots acquired by the missile versus 
range . 
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Figure 13 presents the probability of a missile 
fragment hitting the ship versus the missile destruction 
range. It is clear that the probability increases 
dramatically when the missile is destroyed closer to the 
ship. For a missile, approaching with Mach 2 speed and 7 m 
altitude destroyed by the Phalanx Gun at 200 m from the ship 
and breaking into 20 pieces, there will be about 10 of them 
that finally hit the ship. 




Figure 13. Probability of a single missile fragment to hit 
the ship versus the missile destruction range. 
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b) The Rolling Airframe Missile (RAM) 

Another relatively new weapon system is the 
American/German Rolling Airframe missile (RAM) that was 
first introduced in 1993 to counter the anti-ship cruise 
missile threat. It is a lightweight quick reaction, high 
firepower guided missile weapon system that provides another 
layer of hard kill defense and fills the gap between the 
Phalanx CIWS and the longer-range anti-missile missile 
system as the Nato Sea Sparrow. 

RAM comes with a 21-round guided missile launcher, 
which provides high firepower. The supersonic missile 
demonstrates an excellent kill performance resulting from 
its high maneuverability and accuracy as well as from its 
effective warhead and fusing. Moreover, the latest missile 
upgrade RAM Block 1 Infrared (IRMV) exhibits multiple 
guidance modes, reduced IR and Radar signatures, 
countermeasures capability, and improved performance in 
degraded weather conditions such as heavy rain or fog. The 
currently produced RAM missile (Block 0) has been fired in 
105 flight tests resulting in 100 successful intercepts. 
Its reliability has been sustained in extensive field tests 
under high stress conditions and realistic scenarios thus 
demonstrating that RAM is one of the world's most advanced 
ship self-defense weapon systems. Another system that is 
also under development is the Raytheon's SEA RAM, which 
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combines elements from both Phalanx and RAM CIWS. An 11- 



missile round RAM launcher assembly loaded with RAM B1 
guided missiles replaces Phalanx's 20mm gun. SEA RAM 
inherits RAM's accuracy, extended range and maneuverability 
combined with the Phalanx's Block IB's high resolution 
target detection search and track sensor systems and 
reliable quick-response capability. SEA RAM is an 
affordable capability upgrade and especially attractive 
option for those navies that have already deployed the 
Phalanx, like the Hellenic Navy. It fits the exact 
footprint of the Phalanx, uses the same power and requires a 
minimal ship modification. 

Nevertheless, anti-ship missile technology evolves 
extremely fast, ASMs have not yet reached the peak of their 
evolution. Future high performance threats are expected to 
be stealthier, faster, and lower flying missiles able to 
perform unique maneuvering profiles in order to defeat ship 
defenses. Since a defensive missile is required to have 
substantially superior speed, guidance system, and 
maneuverability than the incoming threat, it is a matter of 
time that RAM or any current missile based inner layer 
defense system will become outdated and ineffective. The 
defensive - offensive missile technology race is clearly won 
by the offensive one. 
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D. DEFENSE AT THE SPEED OF LIGHT 

So far, only kinetic kill methods have been used to 
project power to the target, generally in the form of 
projectiles, rockets and missiles. Kinetic energy weapons 
require a finite flight time to intercept a target. This 
time turns out to be crucial when a threat is approaching 
with supersonic speed, making the reaction time for the ship 
minimum and leaving no chance for reengagement. 




Figure 14. HEL Engagement Scenario. (From Ref. 7). 
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High Energy Lasers (HELs) are now being developed to be 
used as weapons delivering lethal power at the speed of 
light. A high power optical beam traveling with the amazing 
speed of 300,000 km/sec is now available to boost ship's 
self-defense. The HEL provides an effective engagement 
range of approximately 10 km and a practically zero flight 
time out to the target. In the case of a typical short- 
range (10 km) engagement with a modern defensive missile 
system, the time between the employment and threat 
destruction can be as big as 10 seconds. In the 
unfortunate case of a miss there is no time for 
reengagement. On the other hand it takes only 30 psec for 
the laser beam to reach the target and about 2 seconds to 
destroy it. Such a short time provides the benefit of a 
"second chance", if there is a miss and moreover the 
advantage of engaging multiple threats. 

An example of an HEL is MIRACL (Mid- Inf ra-Red Advanced 
Chemical Laser) , which has been developed at White Sands 
Missile Range in New Mexico. MIRACL is a deuterium fluoride 
laser that operates at 3.8 pm wavelength, and has achieved 
an output power of 2.2 MW. At tests conducted in February 
of 1996 under the Nautilus program, MIRACL successfully 
engaged and destroyed a short-range rocket in flight using 
the Sea Lite pointing and tracking system. 
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Undoubtedly, the HEL could enhance substantially a ship's 
self-defense capability, and HEL development should be 
considered seriously to provide an alternative to current 
short-range defense weapons. It is the weapon of the 21 st 
century and will always be effective. Its instantaneous 
response and amazing capabilities can never be outdated, as 
the speed of light will never be exceeded by the performance 
of any future advanced missile. 

E. LASER BEAM LETHALITY 

The HEL weapon concept is fundamentally different from 
that of the kinetic energy weapon. A kinetic energy weapon 
requires a finite time to deliver its warhead or mass to the 
target, but causes damage instantly. On the other hand, the 
laser beam arrives at the target instantly, but it takes a 
finite time to cause damage. The damaging effect of the HEL 
is a thermal process nearly identical to cutting a peace of 
metal with a blowtorch. The question posed is how lethal a 
laser beam has to be in order to shoot down a flying 
missile . 

Although it is preferable, it is not necessary to cause 
the ignition of the warhead in order to destroy a flying 
target. A hole of the dimensions of a child's fist (3cm x 
3cm x 3cm) is generally enough to cause severe damage to the 
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Furthermore, 



missile's electronics. Furthermore, this damage 
significantly alters the missile's aerodynamics resulting to 
instability, overheating and fall. The power required to 
make such a hole depends on the missile's material. As a 
rule of thumb, we can say that it takes a few electron-volts 
(eV) energy to break a chemical bond, and separate atoms. A 
typical 1 liter chunk of material to be removed from the 

missile consists of about 10 25 atoms, so that about 3 MJ 
energy is needed to create such a hole. Limiting the 
delivering time of this energy to 2 to 3 seconds, the needed 
power should be about 1 MW. After considering diffraction, 
the laser beam is focused at the target in an area of 

100 cm so that the intensity is about 10 kW/cm . In order 
to get a feeling of the magnitude of this intensity, we can 

mention that 35 W/ cm on the human body results in third 
degree burns with destruction of the dermis and epidermis. 
It must be noted that delivering time cannot be diminished 
without a limit. Smaller delivery times, which imply very 
high power, may cause detrimental, nonlinear effects such as 
thermal blooming. These effects form a lens in the 
atmosphere and spread the laser beam diminishing the power 
at the target dramatically. Thermal blooming will be 
discussed in the following section. 
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F. 



PROPAGATION OF THE LASER BEAM THROUGH THE ATMOSPHERE 



1 . Transmittance 

The atmospheric components (elements and compounds) 
strongly affect the propagation of the Laser beam and can 
dramatically diminish the required optical intensity at the 
target. There are three primary processes that affect the 
transmission of radiation through the atmosphere: 
absorption, scattering, and random refractive index 
fluctuations (turbulence) . 

Absorption is a quantum process where an atmospheric 
molecule absorbs the energy from some incident photon and 
occurs when certain types of particles interact with the 
propagating electromagnetic radiation. Water, vapor, carbon 
dioxide, ozone, nitrous oxide, and carbon monoxide are the 
major radiation absorbers. Scattering is the redirection of 
radiation by particles in the air. Photons collide with 
these particles, whereas the photon energy is reradiated in 
all directions. Turbulence describes the time-space varying 
temperature and pressure inhomogenuites, which result in 
considerable index of refraction fluctuations. The index of 
refraction variations cause the direction of light 
propagation to bend in various directions. 



34 



Absorption and scattering are grouped together under 
the topic of extinction. Extinction is defined as the 
attenuation in the amount of radiation passing through the 
atmosphere. Absorption and scattering are described 
respectively by the absorption coefficient £ and scattering 
coefficient p, while extinction is described by the 
extinction coefficient p = £ + p [ km ] . Transmittance (Tr) 
through the atmosphere over a distance x is defined as 

— p*x 

Tr = e , and the irradiance at a target at distance x is 

I ( x ) = Io e ^ = Io Tr • Both s and p, and therefore the 
transmittance Tr, strongly depend on the radiation 
wavelength . 

The atmospheric transmittance over a wavelength range 
is extremely crucial to directed energy because it 
determines the choice of the suitable propagating 
wavelength. Transmittance depends on the composition, 
density and pressure of the atmospheric components together 
with the optical line strengths and line widths of all the 
spectroscopic transitions of the molecules and the number, 
size and composition distributions of the particles. 

This complex problem has been described by computer 
modeling codes such as LOWTRAN, FASCODE and MODRAN. MODRAN 
and FASCODE are currently the most widely used and have the 
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capability to handle arbitrary geometries and view 
conditions, while they allow the user to select from several 
model atmospheres, aerosol profiles, and other weather 
conditions. A graph of a typical atmospheric transmission 
in a maritime aerosol model is shown in Figure 15. Note the 
strong dependence of the transmittance on optical 
wavelength. There are specific windows of high 
transmittance where it would be preferable to operate a 
laser. To achieve a wavelength in one of these windows, a 
laser must be tunable. This feature is unique to an FEL; 
other lasers are confined to specific wavelengths. 




Figure 15. Atmospheric Transmittance. (From Ref. 5). 
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2 . Nonlinear Effects 

A high power optical beam traveling through the 
atmosphere is also subject to a "nonlinear" effect where the 
propagation medium is modified by the radiation. This 
nonlinear effect is called thermal blooming, and takes place 
when the absorption of the beam energy by the atmosphere 
causes a local temperature rise followed by a density 
change. Since the refractive index depends on the 

atmosphere density, a local refractive index gradient across 
the laser beam cross-section forms a lens along the beam 
path. This effect can result in spreading, distortion, and 
bending of the laser beam, while it generally reduces the 
power density at the target. In fact, it can limit the 
irradiance to a value independent of the transmitted power. 
Nonlinear effects depend on wavelength, pulse length, 
intensity, and such medium properties as composition and 
density. 

As a result, it is not possible to increase the laser 
power without limit, and continue to have the laser beam 
propagate through the atmosphere. Nonlinear effects can 
lead to dispersion of the beam power, and even complete 
blockage . 
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G. THE FREE ELECTRON LASER (FEL) 

Among the various types of existing lasers, the FEL 
turns out to be an excellent candidate for a ship self- 
defense weapon system. It can be easily integrated in the 
ship's structure and become a major weapon system. The 
electrical power required for the FEL's operation can be 
provided by the ship's power system or by large capacitors. 
The weight and volume of all the FEL components together are 

approximately 150 tons and 140 m 3 respectively [Ref. 8], 
and can be accommodated in a ship the size of a frigate or a 
destroyer. 

The FEL offers significant advantages over other types 
of lasers. The continuous tunability over a wide range of 
wavelengths allows the opportunity to select the frequency 
that best propagates through the atmosphere. FELs do not 
produce hazardous and toxic wastes, which can be life 
threatening for the ship's crew as do chemical lasers. 
Finally the FEL shows the possibility of high wall plug 
efficiency (~10%), making good use of the ship's valuable 
power resources. The fact that it operates with electrical 
power, equips the FEL with an infinite magazine, enabling it 
to minimize reaction time and engage multiple targets. 



38 



IV. SIMULATION RESULTS 



A. THE TAPERED UNDULATOR 

At normal saturation, the electrons lose enough energy 
to drop out of the gain spectrum bandwidth. At resonance we 

have v = 0 and A = A 0 (l + K f / 2y 2 and is only near resonance 
that the coupling between the electron beam and the light 
wave is significant. If A 0 or K are decreased as y 
decreases, as electrons lose energy, resonance would be 
restored, and the interaction could continue to higher 
optical powers. Increasing A 0 or K also extends the FEL 
saturation limit. Decreasing or increasing, tapering or 
inverse tapering, the undulator wavelength A 0 or K along z, 
modifies the simple pendulum equation to include an 
additional torque, or phase acceleration 5, so that 

° OO 

v = £ = 5 + jof| cos(£ + cp) . If 5>0, we have the case of the 
"tapered" undulator, while for 5<0, the undulator is 

"inverse tapered". The value of 5 is given by 

5 = -2nNAA 0 / A 0 when the undulator wavelength is decreased, 

or by 5 = -4nN K 2 AB / B (l + K 2 ) when the undulator field 
strength is decreased. 
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Figure 18. Tapered Undulator, 5-6n. Figure 19. Tapered Undulator, 5=8n. 





Figures 16-19 show the computer-simulated gain 
spectrum, as in Figure 6, of an FEL with dimensionless 
current j=6 and N=41 periods, for various tapering rates 5. 
We can see that as the tapering rate 5 increases, the gain 
spectrum is no longer antisymmetric about resonance, but 
distorts as the bandwidth broadens. The gain in weak 

optical fields ( a 0 « n ) is bigger in the untapered case, 
while in the strong field regime the gain is increased 
because of taper. Figure 20 presents the peak gain of an 
FEL with j=6 as a function of taper, for three different 
optical field strengths a 0 . In weak optical fields, a 0 « 0 

when the magnitude of 5 increases, the peak gain decreases. 
The peak gain curve is symmetric in 5 with peak value at 
5=0. When the optical field becomes stronger, the shape of 
the curve is distorted; the peak value decreases and shifts 
to negative values of 5. Figure 20 shows how the tapered 
undulators have higher gain in strong optical fields than 
the untapered case, 5=0. In strong optical fields a 0 = 40 

with j=6, the maximum gain is 9% achieved with 5=-4n. For 

5=0, the maximum gain in strong fields is only 4%. For 

positive values of 5 and a 0 =40, the maximum gain of 8% is 
achieved with 5=10n. 
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Figure 20. Plot of peak gain versus tapering rate. 

Figure 21 presents a plot of the value of phase 
velocity v 0 giving peak gain in weak fields as a function 

of the tapering rate 5. In weak fields and for moderate 
taper rate (0<5<10n) , the gain peaks at phase velocity 
v 0 =2.4-5/2 . Larger taper rates and strong fields alter 
this relation. 
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Figure 21. Position of Peak Gain versus taper. 



B. MULTIMODE EVOLUTION OF SHORT PULSES 

Multimode behavior is simulated by introducing a short 
parabolic electron pulse shape into the optical resonator 
and following the evolution of the optical pulse and optical 
mode as a function of the number of round trips n, through 
the resonator. Picosecond long electron pulses are often 
used in the FEL oscillator driven by an RF accelerator. At 
resonance, exactly one wavelength of light X passes over 

an electron as the electron passes through one undulator 
wavelength A 0 . Over the whole undulator length, N 
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wavelengths of light pass over the slower electrons. We 
define the "slippage distance" as NX and is the distance 

that the electron pulse slips behind the faster optical 
pulse as they travel through the undulator. It is the 
characteristic length over which electrons and light can 
exchange information during the interaction in the 

undulator. When the electron pulse length is comparable to 
the slippage distance, then short pulse effects are 
important in the FEL interaction. 

As short electron pulses enter the undulator, short 
optical pulses start from spontaneous emission and bounce 
between the resonator mirrors separated by a distance S, 
which is greater than the undulator length L. At t = 0 the 
electron pulses from the RF accelerator enter the undulator, 
while the rebounding optical pulses arrive to overlap the 
electron in time intervals of 2S/c, the round-trip bounce 
time of light in the resonator. We define "desynchronism" 
d = — 2 AS/ (NA) as the displacement between the electron and 
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Figure 22. Short Pulse Evolution with 5=4n and d=0.02. 

Figure 22 shows a sample simulation result after n=2000 
passes with taper 5=4n. Several parameters are printed 

across the top; the peak current j, the resonator loss 
factor Q, the pulse width o z , also 5, and N. The upper 

graphs give the optical field shape |a(z, n)| , the optical 
power spectrum P(v,n), and the electron spectrum f(v,n) at 
the final pass. The shading in the middle graphs show how 
these quantities have evolved with n. The scale ranges from 
0 in black to the maximum value shown as gray. On the 
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bottom left, the longitudinal profile of the current 
density, or electron pulse j(z-x) is shown for reference at 
dimensionless times t=0 (black) and i=l (gray) . The 
calculational window is 10 slippage distances long, and 
electron pulse slips behind one slippage distance in that 
window. The z parameter in this Figure is dimensionless 
( z/NA ) . The bottom center graph shows the weak-field gain 
spectrum for reference, and the right bottom graph shows the 
evolution of the total power P(n). The power is the square 
of the dimensionless optical pulse amplitude a integrated 
over the pulse length. In this particular case, after going 
through some transients, the FEL reaches a steady state 
within approximately n=100 passes. Note the sudden shift of 
the optical power spectrum at n« 100; as n increases, the 
optical field strength grows and distorts the gain spectrum. 
As seen in Figures 16-19, the gain peaks at larger values of 
v as the optical field a becomes large, causing the power 
spectrum to shift to larger values of v. 

Gain in weak fields is described by G(i) = jv 0 T 4 / 12 so 
that there is no gain at the beginning of the interaction. 
This delay in gain is called "lethargy". At later times the 
electron pulse becomes bunched and the light pulse is 
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distorted because gain is preferentially deposited on the 
trailing edge of the pulse. The centroid of the light pulse 
travels slower than c. At exact synchronism (d=0) the 
electron and light pulses enter the undulator 
simultaneously. Surprisingly, the steady-state power of the 
FEL in this case is 0 because the light pulse drifts away 
from the electron pulse over many passes. Figure 23 shows 
the pulse evolution of such a case where the optical power 
evolves to zero after being started with an artificial pulse 
of light. Introducing a small value of desynchronism, we 
compensate for the slower speed of the light pulse and we 
find a better coupling of electron and optical pulses. 
Considerable power is achieved in the steady state, as shown 
in Figure 22. 
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j=10 a z = 1.0 N=41 5=4n Q=10 d=0 




-5 2 5 -47 v 47 0 n 2000 



Figure 23. Short Pulse Evolution with 5=4n and d=0. 



C . TRAPPED-PARTICLE INSTABILITY 

In strong optical fields, electrons can become trapped 
in deep potential wells in phase space. The oscillation of 



the 


bean current 


that is 


trapped in these wells can dr 


ive 


the 


carrier wave 


unstable 


and cause sideband frequencies 


to 


grow 


from noise. 


As the 


FE1 reaches high-power saturati 


on. 



the height of the separatrix is large and many electrons are 
crapped in the closed orbits of phase space. Trapped 
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electrons execute synchrotron oscillation with the 
synchrotron or trapped-particle oscillation frequency 

v s ^(jaj 2 - 5“)^. Sidebands can appear around the fundamental 
at v 0 ± v s , and are shifted away from the fundamental 
wavelength by AA / A = v s / 2nN . The general features of 

trapped particle instability depend only or. j and the loss 
factor Q. Increasing either j or Q can result in the 
trapped particle instability. The trapped particle 
instability usually occurs for small desynchronism, and has 
mixed effects in the FFL. Often, the experiment would like 
a high power laser that has a narrow spectrum and no 
sidebands. 3ut the presence of sidebands means that high 
power has been obtained and the FF1 is otherwise operating 
well. In fact, as more sidebands develop, even more power 
and higher efficiency are attained. Figure 24 presents the 
features of the trapped particle instability. The optical 
pulse (upper left) is short and partially modulated with 
sharp spikes. A small sideband is also evident at the 
optical power spectrum ?(v,n). 
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Figure 24. Trapped Particle Instability. 



D. SIMULATIONS OF THE TJNAF FEL WITH A TAPERED UNDULATOR 

Experiments using the TJNAF FEI have explored operation 
with both the tapered and inversely tapered undulators (Ref. 
15 and 16) . I present here the results of numerical 
simulations using the TJNAF experimental parameters. The 
TJNAF FEL has N=41 undulatcr periods of length A 0 = 2.7 cm 

each. The ur.dulator parameter is K= 0.S3, the resonator 
less factor is Q=10 (10% loss per pass), and the linear 
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taoer 



rates 



are 



AK/K = 5% (5 = 4n) 



. 5 % (5 = 6n) , 



and 10% (5 = 8n ) . Multimode simulations describe the 
evolution of short pulses in the far infrared, and show how 
positive tapering affects single pass gain, steady state 
power, and the electron energy spread as a function of 
desynchronism. 



1 . 34.5 MeV Electron Beam 





n 


this 


case, the FE1 


was ope 


rated with 


short 


0.5 ps 


elect: 


ron 


pul 


ses (length l e 


= 150um ) 


of total < 


energy 


F=34 . 5 


MeV, 


wit 


h a 


0.25% energy 


spread. 


The peak 


curre 


nt was 


1=50 


A, 


and 


the resulting 


optical 


wavelength 


is A 


= 6ym . 



From these values, the dimensionless parameters where 
determined: 1=10 (dimensionless current), and 
a, = l e / NA = 1.0 (electron pulse length). The goal was to 

determine the steady state power, the steady state gain in 
•weak fields, as well as the electron energy spread as a 
function of both 5 and d . 
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a) Steady State Power 



Figure 25 summarizes the results of many 
simulations obtained using numerous values of desynchronism 
from d=0 to d=0.4, for caper rates 5= 0, e4n, -6n, and -8n. 
The number of passes n in each simulation was sufficient for 
the optical power to evolve to steady state. The steady- 
state dimensionless power is plotted versus desynchronism, 
for different tapering rates 5, giving rise no the 
desynchronism curves. 




Figure 25. Desynchronism Curve for various tapering rates. 
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For values of d < 0.004, the laser does not operate. 
3ut at slightly large value of d, the power rises sharply. 
For small d, the gain is small so that a large number of 
passes n is required to achieve steady-state power. For 
5=0, the sharp peak in power at small desynchronism is 
accompanied by the trapped-particle instability. This 
effect is reduced as taper increases, since tapering tends 



to reduce 


the 


side bands. 


As the tapering 


rate is 


increased. 


the 


steady-state 


power decreases and 


the peak 


power of 


each 


curve moves 


slightly to larger 


values of 



desynchronism. The conventional ur.dulator is more efficient 
than the tapered undulators. Furthermore, higher tapering 
rates reduce the range of desynchronism values where the 
laser works. For no tapering, the FEi works for values of 
desynchronism as big as d=0.38. However, a tapering rate of 
5=8n reduces the range to only Ad=0.19. For 5=0, the power 
peaks at d=0.CC4 with the trapped-particle instability. 
This is expected since the sidebands contribute 
significantly to the final power. For 5=8n, the steady- 
state power is 50% smaller and peaks at d=0.015. For larger 
d, the power diminishes significantly as the taper 
increases. The fiat sections of the desynchronism curves 
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since they are 



are desirable regions no operate the FEL, 
predictable and stable. Limit-cycle behavior is observable 
for 5=0 between d=0.C2 and d=0.07, and for 5=8n between 
d=0.025 and d=0.03. For these regions, shown on the graph 
as large circles, only the peak power of the steady state 
power is shown. Limit-cycle behavior causes oscillations 

of the power ?(n), the power spectrum P(v,n), and the 
electron spectrum f(v,n). The steady-state power oscillates 
periodically by as high as 50%, between the peak and the low 
values. Limit-cycle behavior occurs when trapped particles 
in strong fields combine with short optical pulses. The 
modulation caused by the oscillation of the trapped current 
continually modifies the shape of the short optical pulse. 
The different pulse shapes have different powers P and 
spectra P(v), causing oscillations as subpulse structures 
march through the optical pulse envelope. Figure 26 shows 
an example of limit cycle behavior for 6=8n and d=0.C26. 
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j=10 a 2 = 1.0 



N=41 



5=4n Q=10 d=0 . 12 




2000 



P<v,n) 



f (v,n) 




-2 z 2 -50 V 50 -50 V 50 



-2 z 2 -50 v 50 0 n 2000 



Figure 26. Linit-cycle behavior. 

b) Steady State Gain 

In order to determine the steady- state, weak field 
gain as a function of de synchronism, the same simulations 
were used. Instead of plotting power as a function of n in 
the lower right graph of each simulation output, P(n) is 
replaced by gain G(n), as shown in Figure 27. The 
simulation was run only long enough for gain to reach a 
steady-state value. 



oo 










j=10 a z =l . 0 N=41 5=4n Q=oc d=0 . 12 




-5 z 5 -47 V 47 -47 V 47 




-5 z 5-47 V 47 0 n 30 



Figure 27. FEE Short Pulse and Gain Evolution. 

Since we are interested in the weak field gain, 
the optical field a must be kept less than unity by starting 

with an initial optical field a ; * 10 -: . The loss factor is 

Q = 10‘° ~ » , so that the mirror losses (1/Q) would be 
negligible. The power and electron spectra are narrow, 
since the optical field is small. 
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Figure 28. Steady State Gain versus desynchronism 



In Figure 28, the weak-field, steady-state gain 
is plotted as a function of desynchronism d, for the given 



tapering 


rates 


5. 


Clearly the gain curves 


are 


much 


different 


from 


the 


desynchronism curves and 


peak 


for 


moderate 


values 


of d. 


When tapering increases. 


the 


peak 



gain of each curve shifts to smaller values of d. For no 
tapering (5=0), the peak gain is 85%, at d=0.14. When 5=8n, 
the peak gain drops to only 12%, at d=0.03. The most 
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common feature of the desynchronism curve is that the 
operating range, (Ad), decreases as 5 increases. Note that 
when the gain drops to Q=10%, the final power drops to 0 
because the FEL loss per pass exceeds the gain per pass. 

c) Electron Energy Spread 

When the power reaches steady state, we measure 
the full width Av of the final electron spectrum, which is 
the upper right graph of the computer simulation output (see 
Figure 24). From Av we can determine the fractional 
electron energy spread using Ay/y = Av / 4nN where N=41. In 
the TJNAF EEL design, the electron beam exiting the 
undulator is redirected by bending magnets and fed back into 
the accelerator. This recirculation allows energy recovery 
from the electron beam, and significantly increases the FEL 
efficiency. The induced electron energy spread is crucial 
to this process, because beam recirculation is feasible only 
if the full electron energy spread is less than 6%. In 
Figure 29, the fractional electron spread induced by the FEL 
interaction is plotted as a function of d and 5. The curves 
appear much like the power curves in Figure 25 with initial 
sharp peaks of 7%-8% near d=0.01. However for d >0.025, all 
the energy spreads are less than 6%. 
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Figure 29. Electron Energy Spread. 



Between d=0 and d=0.05, where we have the maximum 
output power, the induced energy spread is significantly 
smaller for the conventional undulator than the tapered one. 
For d=0.025, the energy spread was 4.5% for the untapered 
undulator and around 6% for all the other tapers. The 
tapered undulator does not appear to be desirable for 
recirculating the electron beam. After d=0.05, the energy 
spread drops to less than 4% for all the tapering rates. 
With short optical pulse, the head and tail of the pulse 
have weak optical fields, so that tapering is not optimum 
along the pulse. This is the reason the energy spread is 
greater in the tapered undulators. 
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2 . 



47.5 MeV Electron Beam 



a) 


Steady- State Power 




In 


this 


case, the 


dimensionless 


parameters are 


CJ. 

II 


and 


o 2 = 1.8 . 


The optical 


wavelength was 



A = 3pm. The desynchronism curves were again determined as 

shown in Figure 30. The curves are very similar with those 
using the 32.5 MeV (Fig. 25). However, the operating ranges 
and the final power are smaller because of the smaller value 
of j. For taper 5=8n, the laser barely works. 




0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 

desynchronism 



Figure 30. Desynchronism Curve for 47.5 Mev electron 
beam and various tapering rates. 
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The operating width of the desynchronism curves, 
Ad, as a function of taper has been measured in experiments 
conducted with the IR Demo PEL at Jefferson Lab [Ref 16] . 
The FEL parameters used in the simulations were the same as 
those of the experiments (A = 3pm, N=41, Q=10, j=7 and 

o 2 = 1.8) . In Figure 31, the relative width of the 
desynchronism curve is plotted as a function of taper, using 
both experimental and computer simulation data. The 
relative width is the maximum value of desynchronism where 
the FEL operates for each taper, divided by the width of the 
desynchronism curve of the untapered undulator (5=0) . Both 
simulations and experiments show that the operating range 
decreases as the taper increases. There is good agreement 
between theory and experiment. 




Figure 31. Relative width of desyncronism curve versus taper 5. 
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The trapped-part icle instability and limit-cycle 
behavior were again characteristic features of the simulation 
results. Figure 32 presents the simulation output with j=7 
and c_ = 1.8 (47.5 MeV Energy), 5=0, and d=0.C25. The 
simulation was run only for n=3G0 passes in order to better 
present the details. The oscillation of the total power, 
the optical pulse distortion and the sideband in the optical 
power spectrum are clearly shown. 

j=7 a 2 =l . 8 N=41 5=0 Q=10 d=0 . 025 



lCC(z,n) j P (V/n) f (v,n) 




-5 z 5 -47 V 47 -47 V 47 




-5 z 5 -47 V 47 0 n 300 



Figure 32. Simulation Output for d=0.025 and 5=0. 
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Limit-cycle behavior in the final power was 
observed for 5 = 0, 4n, 8n and not for 5 = 6n. Figures 
33-35 show the desyncronism curves of Figure 30 in order to 
more clearly presentthe details of limit-cycle behavior. 




Figure 33. Desynchronism Curve for 5=0. 
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Figure 34. Desynchronism Curve for S=4n. 

The upper right graphs of Figures 33 and 34 show 



blown up section of the desynchronism curve section where 
the limit-cycle behavior occurs. Since there are periodic 
oscillations of the total power, the maximum and minimum 
power is indicated in the graph. The mean power is the 
average of the maximum and the minimum. The desynchronism 
curves in Figure 30 plotted the mean power. Figure 35 shows 
the desynchronism curve for 5=8n, emphasizing the limit- 
cycle behavior. 
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Figure 35. Desynchronism Curve for 5=8n. 



b) Steady-State Gain 

The same procedure as the one followed in 
paragraph lb was followed in order to determine the steady- 
state, weak-field gain as a function of taper and 
desynchronism. The results are presented in Figure 36 and 
are similar to those with the 32.5 Mev beam energy 
(Figure 29) . Still the operating ranges, as well as the 
gains, are smaller. 
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N=41, 1.8 




Figure 36. Steady State Gain versus desynchronism for 47.5 
Mev electron beam. 



c) Electron Energy Spread 

Figure 37 shows the fractional electron spread as 
a function of d and 5. The results are again similar to 
those with 32.5 MeV electron beam energy. The curves have 
the same trend as the power curves (Fig. 30) with a peak 
near d=0.01. The peak fractional energy spread is 
approximately 1% smaller than for the 34.5 MeV beam because 
of the smaller value of j . 
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7 % n 




Figure 37. Electron Energy Spread for 42.5 MeV Electron Beam. 

Initial peaks of 7% induced energy spread are 
evident at the start up region (d=0.007). For d > 0.015, 
the energy spread drops to less than 6%, implying the 
possibility of beam recirculation. In the region of high 
steady-state power (0. 01< d <0.04) the conventional undulator 
(5=0) has almost 1% less energy spread than the tapered 
(5=4n, 6n) . Thus, the untapered undulator is a better 
candidate for beam recirculation. The short optical pulses, 
with weak fields at the head and tail, cannot be easily 
designed to reduce energy spread. 



67 



E. CONCLUSIONS 

Experiments were conducted with the TJNAF IR demo FEL 
(Ref. 16) in order to test its operation with various taper 
rates 5, and desynchronism values using the following 



parameters : 



Undulator 

Periods 


N 


41 


41 


Wavelength 


A 


3 pm 


6 pm 


Electron Beam 
Energy 


E 


47.5 MeV 


34.5 MeV 



Table 3. Parameters used in the FEL experiments. 



In this thesis, computer simulations were used in order 
to explore the operation of the TJNAF FEL with the 
conventional and several tapered undulators using the same 
parameters as in experiments. The corresponding 
dimensionless parameters used in the simulations are listed 



in Table 4 . 



Electron Beam 
Energy 


E 


47.5 MeV 


34.5 MeV 


Dimensionless 

Current 


j 


7 


10 


Resonator 
loss factor 


Q 


10 


10 


Pulse length 


O z 


1.8 


1.0 



Table 4. Parameters used in the simulations. 



The desynchronism curves were determined at each 
electron beam energy, and present a sharp peak at small 
values of desynchronism accompanied by the trapped-particle 
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instability. Trapped-particle instability was reduced for 
larger values of taper, since tapering suppresses the 
sidebands. Limit cycle behavior was also observed for most 
tapering rates. Although the general shape of the curves is 
very similar for each beam energy, the final power and gain 
were significantly smaller for the 47.5 MeV electron beam 
than for the 34 . 5 MeV electron beam, because of the lower 
value of j . 

The tapered undulator turned out to be less efficient 
than the conventional undulator, demonstrating reduced gain 
and power, and increased energy spread to the electron beam. 

The width of the desynchronism curve decreased as the 
tapering rate is increased, in agreement with the 
experiments. A publication is being prepared comparing the 
experimental (Ref. 16), and simulation results. 
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V. EXPERIMENTAL RESULTS OF DAMAGE INDUCED TO MATERIALS BY 

THE TJNAF FEL 

A. LASER MATTER INTERACTION 

The interaction between lasers and matter is a 
complicated issue generally described by non-linear and 
unpredictable effects. The laser beam has the unique 
ability to deliver very high power per unit area. When high 
power laser radiation falls on a target, the part of the 
beam that is absorbed begins to heat the target surface very 
rapidly to its melting temperature. This melting process 
then penetrates progressively into the material. Many 
physical processes govern the damage caused to the material 
by the laser power including power absorption, power 
reflection, heat conduction, and heat diffusion. 
Furthermore, a large number of parameters play a major role 
in these processes such as material density and heat 
capacity, as well as the irradiation wavelength, power 
density, peak power, and possibly pulse characteristics. 

A good knowledge of these mechanisms helps one 
understand the capabilities and limitations of the directed 
energy, and allows better control of the damage caused by 
the laser. Controlled damage has many industrial 
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applications such as the creation of thin coatings, 
electronic component fabrication, precise drilling, cutting, 
etc. However, when using the high-power laser as a weapon 
to shoot down missiles, precision and symmetry of the 
damaged area are not the issue. Our goal is to cause the 
maximum possible damage on a rapidly moving missile; this 
requires precise tracking of the laser beam through a 
turbulent atmosphere 



B . SCALING 

As stated in chapter III, an intensity of 10 kW / cm 2 
over a 6 cm radius spot for a few seconds dwell time is 
needed to burn an adequate hole in a missile. In these 
damage experiments, the TJNAF FEL used only several hundred 
watts. In order to achieve the desired intensity of 10 

kW / cm , the laser beam has to be focused with a lens to a 
much smaller spot size. For average power P and spot area 
A, the resultant intensity is I = P / A . Given a specific 
power, we can adjust the spot size to achieve the desired 
intensity. Studying a spot radius on the order of 0.1 cm 

and using 500 W power, we can achieve 10 kW / cm intensity. 
The goal of this scaling is to develop guidelines that will 
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reliably predict the damage caused by a high power laser 
over a large area. 

Nevertheless, the spot size cannot be arbitrarily 
small, or the scaling will not work. Each material is 
characterized by a parameter called the thermal diffusion 
length D, which represents the distance required for the 
temperature to drop by a factor of 1/e, and determines the 
ability of the material to absorb and transport heat. If D 
is greater than the laser beam spot size diameter, the heat 
deposited by the laser beam will diffuse away in less time 
than it takes the material to melt. In order to obtain 
effective heating, and cause melting of the material 
irradiated, the thermal diffusion length should be smaller 
than the beam diameter. Also, the damage hole cannot be 
much greater than the laser beam radius. Material cannot 
easily escape a deeply damaged hole, and the damage 
mechanism may be complicated by the material ejecting the 
hole and flowing through the beam. 



C . EXPERIMENTAL PROCEDURE 

The sample materials used in the experiments were 
Polyimide Fiberglass and F2 Epoxy. The samples provided by 
the Naval Research Laboratory are the same as those used in 
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similar experiments conducted in March 1999. The results of 
the March 1999 irradiations, which involved three holes in 
each sample labeled 1,2, and 3, were presented in Reference 
[14], and some of them will be mentioned in this thesis. 
The results presented by this thesis were obtained by two 
different experiments conducted in August 1999 and March 
2000. The results from both experiments will be presented 
together in order to make the analysis more comprehensive. 

1. August 1999 Experiment 

Each sample was irradiated 12 times (4 sets of 3 
irradiations with the same parameters) with a FEL beam of 
wavelength A=3.1 pm, pulse repetition frequency 18.7 MHz and 
average power P=100W±5W. The Polyimide Fiberglass was 
irradiated first, followed by the F2-Epoxy. Measurements 
were made with samples placed downstream of a calcium 
fluoride lens with a measured back focal length of 137.6 mm 
for 3 pm wavelength. A camera was set up to observe the 
front and back surface of the samples. Two sets of three 
irradiations were first made. The average intensity was 500 

W / cm 2 , which was achieved by focusing the beam to a spot 
of 0.25 cm radius. Three identical irradiations were made 
initially with no airflow. Then three more irradiations were 
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made with a wind speed of 60 mph across the front face of 
the samples. Then, adjusting the beam radius to 0.087 cm, 

the intensity was set to 1=10 kW / cm 2 and the same set of 
measurements were repeated. 

2 . March 2000 Experiment 

The goal of these measurements was to maintain the 

average intensity of 10 kW / cm 2 , but using higher laser 
power with a larger spot size. Each sample was again 
irradiated 6 times (2 sets of 3 irradiations) with a FEL 
beam of wavelength X=3 . 1 pm, pulse repetition frequency 
37.425 MHz and average power P=500W±10W. The Polyimide 
Fiberglass was irradiated first, followed by the F2-Epoxy. 
Measurements were made with samples placed downstream of a 
calcium fluoride lens with a measured focal length of 235.7 
mm for 3 pm wavelength. A camera was again set up to 
observe the front and back surface of the samples. 

Irradiations were made with 1=10 kW / cm average intensity, 
which was achieved by focusing the beam to a spot of 0.13 cm 
radius. Three identical irradiations were made with no 
airflow, and then three more irradiations with a wind speed 
of 85 mph across the front face of the samples. 
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D . RESULTS 

1 . Sample #1 - Polyimide Fiberglass 

The sample had dimensions 11.4 cm by 10.1 cm with 2 mm 
thickness. In Figure 38, we see a photo of the front side 
of the Polyimide Fiberglass sample after all sets of 
irradiations. Irradiations labeled by numbers 1,2 and 3 
were conducted in March 1999. Irradiations labeled 4 to 15 
were conducted in August 1999, and 16 to 21 were conducted 
in March 2000. All irradiations were done three times with 
the same parameters in order to get more accurate 
measurements. The actual results came from the mean value 
of the three measurements. Irradiations 7,8,9,10,11,12,16 
17 and 18 were done with no airflow, while in 
4,5,6,13,14,15,19,20, and 21 there was wind present. In 
Figure 39, we see a photo of the backside of the Polyimide 
Fiberglass sample showing that all of the irradiations with 

intensity of 10 kW / cm completely penetrated the sample. 
On the other hand, the irradiations with an intensity of 

500 W / cm 2 (irradiations 10 to 15 of Figure 38) did not 
penetrate the sample. Irradiation results are summarized in 
Table 5. Each row represents a set of three irradiations. 
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The hole diameters have variation of 15% in each set of 



three holes. This is due to the slight variation of the 
exposure time, which is on the order of 1-2 seconds. 

The presence of airflow resulted in a 15% bigger hole 
diameter and penetration depth rate, than those attained 
without airflow. The irradiations conducted with 500 W 
average power (runs 16 to 21), caused 3.5 times higher 
penetration depth rate than the ones conducted with 100 W, 

and the same intensity ( kW / cm 2 ) . This is probably because 
the higher power allowed larger spot, so that damaged 
material could more easily escape the larger damage hole. 

The damage pattern of runs 16 to 21 is not circular as 
expected but somewhat elliptical, which is very clear in 
Figures 44, 45,46 and 47. The following figures present a 
closer caption of the damage caused by the irradiations. 
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Figure 40. Poiyirr.ide hole 12. 

2 

Intensity 5C0 W / ex , no air. 
Average Power 100 W. 

Spot Radius 2.5 mm. 



*1 ^ 
air. 

Average Power 100 W. 

Spot Radius 2.5 mi. 



Figure 41. Poiyirr.ide hole 

2 

Intensity 500 W / cm. with 




Figure 42. Polvinide hole 


£ . 


Figure 43. Polyimiae hole 5. 


2 

_ntensitv 1C kW / cm , no 


air. 


2 

Intensity 10 kW / cm , with air 


Average Power 100 a. 




Average Power ICO W. 


Spot Radius 0.87 mm. 




Spot Radius C.87 mm. 



8C 








Figure 44. Polyimide hole 18. 

2 

Intensity 10 kW / cm , no air. 
Average Power 500 W. 

Scot Radius 1.25 irn. 




Figure 45. Polyimide hole 20. 

2 

Intensity 10 kW / cm , with air. 
Average Power 500 W. 

Scot Radius 1.25 mm. 
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In 


all cases a 


raised 


lip of melted 


material is 


observed 


around the 


face of 


the entrance 


hole. The 



dimensions of the lip are approximately 0.2 mm height and 
1 mm width. However, when airflow is present the lip tends 
to be smaller, possibly because the airflow removes the 
debris, and the evaporated material that actually 
contributes to the formation of the lip. That would explain 
the bigger diameter of the holes made in the presence of 
wind, since the diameter is measured from the inside part of 
the lip. 

The charred region extends approximately 2 mm around 
the lip when there is no wind. With wind this area is much 
smaller. The charred region around the exit hole extends to 
1 to 1.5 mm in all cases, as the backside of the sample was 
not exposed to the airflow. As exposure time increased the 
radial extent of the damaged area increased and more melted 
material was deposited around the hole. After the 
irradiation stopped, there was a period of almost 3 sec that 
the material is still hot and melting. The wind tends to 
cool down the material, decreasing this time almost in half 
and resulting in less melted material. 

Investigation with a microscope reveals that there is 
no evidence of melted or rehardened material inside the 
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holes . 



It is also evident that the damage is most 



significant in the center of the hole, diminishing radially 
outward, which supports expectation of a Gaussian shape of 
the laser beam intensity. 

The following table presents the results of the 
experiment conducted in March 1999 (Holes 1,2,3). 



Average 

Power 


Average 

Intensity 


PRF 


A 


Spot 

radius 


Wind 


Penetration 

Depth 


(W) 


(kW/cm 2 ) 


(MHz) 


(pm) 


(mm) 


(mph) 


Rate 

(mm/sec) 


100 


10 


37.4 


4.825 


0.87 


NO 


0.9 



Table 6. Irradiation results of March 1999 experiment on 
Polyimide Fiberglass. 

The above results have the same irradiation parameters 
with those of runs 7,8,9 of Table 5 except for the higher 
Pulse Repetition Frequency (PRF) of 37.4 MHz and the longer 
wavelength A=4 . 8 pm. It appears that the shorter wavelength 
in Table 5 (A=3.1 pm) combined with the lower PRF is much 
more effective, resulting in 60% higher penetration depth 
rate than in table 6. The lower PRF apparently causes more 
damage, since it results in more energy per pulse for the 
same average power. 
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2 . 



Sample #2 - F2 Epoxy 



The sample had dimensions 11.5 cm by 10 cm with 1.5 mm 
thickness, attached to a 1.6 cm thick polyethane foam 
backing which is clearly seen in Figure 49. In Figure 48, 
we see a photo of the front of the sample after all sets of 
irradiations. Irradiations labeled 1, 2 and 3 were those 

conducted in March 1999. Irradiations 4 to 15 were 
conducted in August 1999, and 16 to 21 were the last ones 
conducted in March 2000. Following the same procedure as 
with polyimide, all the irradiations were done three times 
with the same parameters in order to get more accurate 
measurements. The results presented here came from the mean 
value of the three measurements. Irradiations 7,8,9,10, 
11,12,16,17, and 18 (Fig. 48) were done with no airflow 
while in 4,5,6,13,14,15,19,20,21 there was wind present. In 
Figure 50, we see a photo of the backside of the sample 

showing that all of the 10 kW / cm irradiations completely 

penetrated the sample. On the other hand the 500 W / cm 2 
intensity (runs 4 to 9) did not penetrate, but caused more 
extensive surface damage due to the bigger spot radius. 
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Figure 48. F2 Epoxy (front view) 




Figure 49. 



F2 Epoxy (side view) . 



85 




The preser.ee of the foam layer at the 
ample made the measurements of the exit 
r.reiiable. A closer caption of the damage 
icures 51-56. 



backside of th 
holes diamete 
is presented i 



Irradiation results are summarized in Table 7. Sac 
cw of the table represents a set of three irradiations 
cnducted with the same characteristics. 
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Figure 53. Epoxy, hole 16. 

Intensity 10 kW/cn 2 , no air. 

Average Power 500 W. Spot Radius 1.25 ran. 




Figure 54. Epoxy, hole 20 

Intensity 10 kW / era", with 
Average Power 5CC K . Spot 



airflow 

Radius 



25 ran. 
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Figure 55. Epoxy, hole 12. 

_r.ter.sity 10 kW / err , r.o air. 

Aver ace Power ICC W. Soot Radius 0.87 nra. 




Figure 56. Epoxy, hole 15 

Intensity 1C kV. / c rr. 2 , with 
Average Power ICO W. Spot 



airflow. 
Radius 0.87 



itra. 
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The hole diameters have variation of 15% in each set of 



three holes. This is again due to the slight variance of 
the exposure time of the irradiations, which is on the order 
of 1 to 2 seconds. The presence of airflow did not seem to 
increase the extent of the damage. The entry hole diameter 
was actually decreased by 10% to 30%. 

The average power of 500 W (runs 16 to 21), resulted in 
2.5 times higher penetration rate than the 100 W power (runs 
10 to 15) , and caused a slightly elliptical damage pattern. 
During the irradiation, it was observed that flames, smoke, 
and debris covered the entrance hole. After the irradiation 
stopped the material was still burning for almost 3 seconds, 
which caused extra charring and melting of the sample. When 
airflow was applied, the time decreased by half. The 
charred region extends approximately 0.5 mm around the 
entrance hole with wind present and 1 mm without wind. 
Examination of the holes with a microscope revealed more 
roughness than the polyimide. This was probably caused by 
deposited debris and charred material. The Gaussian beam 
caused the same damage pattern with the polyimide sample, 
being more intensive at the center of the spot. However the 
penetration rates were observed to be 2 to 3 times smaller. 
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Table 8 presents the results of the experiment 



conducted in March 1999 for the Epoxy sample (Holes 1,2,3). 



Average 

Power 


Average 

Intensity 


PRF 


A 


Spot 

radi 


Wind 


Penetration 

Depth 


(W) 


(kW/cm 2 


(MHz) 


(pm) 


us 

(mm) 


(mph) 


Rate 

(mm/sec) 


100 


10 


37.4 


4.825 


0.87 


NO 


0.1 



Table 8. Irradiation results of March 1999 experiment on Epoxy 
Comparing these results with runs 10 to 12 of Table 7, 
it is clear again that the A=3.1 pm wavelength combined with 
the lower Pulse Repetition Frequency was more effective, and 
resulted in 6 times higher penetration depth rate. 



E. CONLU S I ON S - SUGGESTIONS 

When higher laser power is available in the future, the 
same irradiations should be conducted with even larger spot 
radii in order to compare the results and establish scaling 
laws. It is suggested that thicker samples should also be 
tested in order to determine if the penetration rate stays 
the same after the laser beam has penetrated the material a 
few mm. 
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Airflow did not have a significant effect, possibly 
because of its relatively low speed. It is suggested that 
in future experiments the wind speed should be much higher 
in order to simulate more realistically the conditions of a 
missile flying with 1 to 2 Mach speed (1000 mph) . 

Furthermore, a way to make the measurements of burn 
through times more accurate should be established, as they 
are important in determining penetration depth rates. 

From the analysis of the results, it is apparent that 
the penetration depth rate is proportional to the laser 
intensity. When the intensity is changed by a factor of 20 

(from 500 W / cm to 10 kW / cm ) the penetration depth rate 
is also changed by almost the same factor in both samples. 

The A=3.1 pm wavelength appears to be more effective 
than the A=4.875 pm wavelength. In addition, the lower PRF 
(18.7 MHz) seems to contribute to higher penetration rates 
due to the higher energy per pulse. It is recommended in 
future experiments to keep the pulse energy constant, but 
increase the repetition rate, to see how the burn-through 
rate is affected. It would be useful also to vary the 
wavelength, while keeping the PRF and the energy per pulse 
constant . 
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